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Effects  of  Inhomogeneous  Filtering  and  Buoyant  Forcing 
on  the  Dynamic  Subgrid  Scale  Model 
M.  Bohnert  and  J.  H.  Ferziger 
Dept,  of  Mech.  Eng.,  Thennosciences  Division 
Stamford  University,  Stanford,  California 

Abstract 

The  d3rnamic  subgiid  scale  model  is  used  in  large-eddy  simulation  of  the  stratified  Ekman 
layer.  Results  from  direct  numerical  simulations  (DNS)  for  the  neutrad  Ekman  layer  were 
used  to  test  model  efficacy.  Three  dimensional  test  filtering  was  employed  in  the  dynamic 
calculation.  A  significant  increase  in  accuracy  was  observed  in  the  neutrad  layer  when  this 
filtering  wais  employed,  as  compared  to  filtering  only  in  planes  parallel  to  the  surface.  The 
model  shows  the  ability  to  model  both  unstable  and  stable  regions,  and  its  effect  vanishes 
where  stratification  strongly  quenches  the  turbulence. 


Introduction 


In  large-eddy  simulations  of  turbulence  (LES),  the  largest  scales  are  directly  computed 
while  the  subgrid  scales  (SGS)  are  modelled.  The  latter  are  the  scales  that  cannot  be 
resolved.  The  majority  of  SGS  modelling  closures  make  use  of  eddy-viscosity  parameteri- 
zations;  the  most  popular  is  due  to  Smagorinsky^  in  which  the  eddy  viscosity  is  represented 


UT  =  {CsAf\S\  (1) 

where  A  is  the  filter  length  scale,  C,  is  a  specified  constant,  and  t-Si  =  (25ySy)^^^  is  the 
magnitude  of  the  resolved  strain  rate  tensor. 

While  the  Smagorinsky  model  has  enjoyed  success,  its  application  to  a  wide  variety  of  flows 
has  been  limited  by  the  following  factors: 

(a)  The  viscosity  does  not  vanish  in  laminar  flow, 

(b)  The  model  is  absolutely  dissipative;  that  is,  backscatter  from  small  to  large  scales 
cannot  be  represented, 

(c)  Ad  hoc  modifications  are  often  required  in  regions  of  large  gradients  (e.g.  near  a 
boundary), 

(d)  There  is  widespread  debate  among  researchers  about  the  value  of  C«,  with  an  apparent 
lack  of  a  ‘universal’  value  that  can  be  used  for  a  wide  range  of  flows. 

In  addition  to  the  above  limitations  present  in  engineering  flows,  the  meteorological 
community  has  found  that  the  Smagorinsky  model  is  poorly  suited  to  flows  affected  by 
unstable  and  stable  stratifleation.  The  SGS  parameterization  in  these  flows  commonly 
takes  one  of  two  forms^: 

(e)  Modifications  to  the  basic  Smagorinsky  concept  are  employed;  for  example,  Cs  becomes 
a  function  of  the  local  Richardson  number, 

(f )  The  eddy  viscosity  is  determined  by  the  local  value  of  the  SGS  turbulent  kinetic  energy, 
65.  This  reqriires  a  solution  of  am  additional  equation  for  the  transport  of  £5. 
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Recently  Germano  ei  al.^  proposed  a  dynamic  subgrid  scale  model  which  calculates 
the  Smagorinsky  coefficient  Cs  interactively  as  a  function  of  both  space  and  time.  The 
coefficient  is  computed  via  the  application  of  a  test  filter  with  length  scale  A  >  A  to  the 
LES  field.  The  resolved  information  contained  between  these  two  length  scales  is  used  to 
calculate  C5. 

The  dynamic  model  has  been  successful  in  both  incompressible  ztnd  compressible  homo¬ 
geneous  flows,  as  well  as  incompressible  turbulent  and  trzmsitional  channel  flows^’^.  None 
of  the  limitations  (a)-(c)  listed  above  are  evident  (limitation  (d)  is  no  longer  meaningful 
since  Ca  is  not  constant).  In  turbulent  channel  flow,  however,  the  test  filter  has  only 
been  applied  in  two  dimensions;  these  are  the  planes  parallel  to  the  wzdls  on  which  the 
computational  grid  is  uniform. 

The  purpose  of  this  study  is  to  investigate  the  applicability  of  the  dynanoic  SGS  model  in 
the  LES  of  a  stratified  external  wall-bounded  flow,  the  turbulent  Ekman  layer.  The  efficacy 
of  three-dimensional  test  Altering,  including  an  inhomogeneous  direction  in  which  the  grid 
is  not  uniform,  is  examined  and  compared  to  the  previously  used  two-dimensional  approach. 
A  temperature  profile  with  unstable  and  stable  re^ons  will  be  imposed  on  the  flow  to 
introduce  buoyantly-forced  turbulent  dynamics.  This  stratified  flow  is  an  idealization  of 
a  meterological  application,  the  Convective  Boundary  Layer  (CBL),  and  will  be  used  to 
examine  the  applicability  of  the  dynamic  model  in  flows  where  parameterizations  (e)  and 
(h)  above  had  previously  been  employed. 

Problem  Formulation  and  Numerical  Method 


We  consider  the  flow  of  a  stratified  viscous  fluid,  driven  by  a  uniform  pressure  gradient 
(geostrophic  conditions  aloft)  and  experiencing  steady  system  rotation.  Density  variations 
are  assumed  to  have  no  dynamic  influence  other  than  to  provide  a  temperature  dependent 
body  force  (the  ‘Boussinesq’  approximation).  A  three  dimensional  filter  that  commutes 
with  differentiation  is  applied  to  the  governing  equations;  filtering  the  nonlinear  terms 
results  in  a  subgiid  scale  stress  and  heat  flux,  rg  and  respectively  : 


at 


1  dp 


d'^uj  drs^ij 


— (U.U,)  =  g- 


(3) 


3 


=0 


d9  d  ei^e  dqs,j 

^  dxj 

Here,  fJ,-  is  the  applied  rotation  and  Gj  is  the  free  stream  velocity,  which  edong  with  the 
acceleration  of  gravity,  g,  define  the  imposed  uniform  pressure  gradient 


-Q^  —  — 2poo€jjifljGj|.  —  Poog 


The  pressure  p  in  (5)  is  thus  the  deviation  &om  P  defined  by  (6).  We  have  adopted  the 
meteorological  convention  of  using  (i,y,r)  =  (*i,X2,X3)  as  the  axes  in  the  streamwise, 
spanwise,  and  wall-normal  directions,  respectively.  The  nondimensional  temperatiire,  B,  is 
defined  as 

r  —  Too 


where  T  amd  Too  are  the  local  and  reference  temperatures,  and  <  ...  >  denotes  averaging 
over  (x,  y)  planes.  In  large-scale  meteorological  applications,  the  ‘potential’  temperature 
is  used  to  distinguish  between  adiabatic  and  diabatic  temperature  change.  At  the  scales 
considered  in  this  idealization,  the  hydrostatic  pressure  variation  across  the  depth  of  layer 
is  small,  and  there  is  a  negligible  difference  between  the  actual  temperature  and  ‘potential’ 
temperature. 

The  SGS  terms  rg  and  qs  in  (3)  and  (5)  are  defined  as 


T’s.ii  = 


=  uiB  —  uiB 


and  will  be  modelled.  The  models  are  the  principal  issue  of  this  paper  and  are  discussed 
below.  The  boundary  conditions  used  are 

Uj  =  0  at  z  =  0  (9.1) 

Ui—*Gi  as  z  —*  oo  (9.2) 

B  =  9q  ot  z  =  0  (9.3) 

^  ►  0  as  z  -+  oo  (9.4) 
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All  dependent  variables  are  assumed  to  be  periodic  in  the  planes  of  homogeneity  (x,y). 

The  solution  technique  employed  is  ftiUy  spectral  in  all  three  Cartesian  dimensions  and 
second  order  accurate  in  time.  The  explanation  of  the  method  requires  some  length,  and  the 
reader  is  referred  to  Coleman,  et  al.^  for  a  complete  discussion.  In  short,  the  semi-infimte 
domain  0  <  z  <  oo  is  mapped  into  a  finite  domain  1  >  »7  ^  0  via  the  tramsformation 
Tj  =  where  Z  is  a  mapping  length  scale.  Fourier  methods  are  used  to  represent 

variables  in  the  homogeneous  (x,y)  directions.  A  family  of  orthogonal  Jacobi  polynomials 
is  used  to  generate  basis  functions  in  the  variable  q.  Therefore,  each  dependent  vzuiable  is 
a  generalized  Fourier  series;  for  example, 

M 

u(x, », ,)  =  YIT,  E  (12) 

i,  jb,  m=Q 

where  each  Hm  is  a  basis  function  unique  to  u  that  satisfies  the  boundary  conditions.  The 
basis  functions  used  for  u,  v,  and  w  satisfy  the  continuity  equation  (4)  everywhere  so  a 
Poission  solution  for  pressure  is  not  required. 

In  Fourier  space,  a  filter  may  be  defined  by  a  wavenumber  ‘cutoff’;  that  is,  wavenumbers 
higher  than  some  value  kc  are  discarded  and  a  finite  series  is  formed.  The  cutoff  filter  has 
a  length  scale  defined  as  A  =  7r/|fcc|)  a^d  it  can  easily  be  shown  that  this  filter  commutes 
with  differentiation.  For  variables  represented  by  polynomial  basis  functions,  the  use  of 
a  finite  series  implies  an  effective  ‘cutoff’  filter  that,  in  general,  does  not  commute  with 
differentiation.  In  this  study,  a  filter  is  defined  explicitly  as  follows 

M 

«(*, ».  =  EE  E  /m  =  (13) 

Is,  Jby  m=0 

where  7  is  a  constant.  We  make  the  assumption  that  this  filter  commutes  with  tj  differen¬ 
tiation  (or  equivalently,  with  z  differentiation);  this  is  justified  by  two  observations.  First, 
the  am  coefficients  decay  rapidly  at  high  m.  Since  any  error  in  a  ‘cutoff’  series  is  directly 
related  to  the  magnitude  of  the  first  few  eliminated  coefficients,  this  error  will  be  small. 
Second,  the  additional  application  of  the  explicit  filtering  factor  /m  reduces  the  importance 
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of  these  coefficients  even  further;  in  the  specific  cases  presented  in  this  paper,  7  fa  0.01, 
which  results  in  the  last  coefficient  of  a  t3rpical  series  (i.e.  M  =  20  )  being  multiplied  by 
fm  ^  0.02. 

At  first,  it  may  seem  simpler  to  apply  an  m  ‘cutoff’  filter  to  the  above  series  rather 
than  use  the  fm  approach.  However,  the  continuity  equation  puts  constraints  on  the 
relationships  among  the  am  coefficients;  a  large,  computationally  expensive  linear  set  of 
equations  must  be  solved  if  M  is  to  be  reduced  and  the  boundary  conditions  satisfied. 
Application  of  the  filtering  factor  fm  retains  both  the  continuity  constraint  and  boundary 
conditions,  and  is  computationally  efficient. 

The  definition  of  this  explicit  filter  in  the  inhomogeneous,  or  z  direction,  allows  one  to 
examine  its  effect  on  dynamic  subgrid  scale  model.  This  will  be  discussed  in  detail  in  the 
next  section. 

LES  Modelling 

The  SGS  terms  defined  in  (8)  are  modelled  as  follows 

rs,ii  -  f  TS,U  =  -2^T5i;  =  CA^|5|?«  (w-D 


17*  96  96 

~  ~  PrT9xj  ~  ~'^9xj 


(14.2) 


The  composite  length  scale  A  is  defined  as  (AiAyA*)^/^.  The  coefficient  C  and  turbulent 
Prandtl  number  Pry  are  computed  by  the  dynamic  procedure  proposed  by  Germano  et  al.. 
In  this  procedure  a  test  filter  (denoted  by  a  tilde,  i.e.  u)  with  length  scale  A  >  A  is 
applied  to  the  LES  flow  fleld.  The  resolvable  stress  and  heat  flux  terms  are  computed  from 

Lij  =  UiUj  —  (15.1) 

(15.2) 

Lij  and  Kj  can  then  be  expressed  in  terms  of  modelled  quantities 


Lij  =  -2C(A^|5|S(y  -  A^|5|Siy)  =  -2CMij 


(16.1) 


(16.2) 


^  (A-|5|.^-aV|^)=  ^ 


Prj> 


'dx. 


Values  of  C  and  Prj*  C2in  be  determined  by  contracting  both  sides  of  (16.1)  emd  (16.2) 
witb  Mij  and  Nj,  respectively,  in  a  least  squares  approach  suggested  by  Lilly^. 

We  have  not  yet  defined  the  test  filter  or  the  spatial  dependence  of  C  and  Pry.  To 
address  this  issue,  we  note  that  past  researchers  who  used  the  dynamic  model  in  flows 
with  an  inhomogeneous  direction,  e.g.  the  wall  normal  direction  in  ttirbulent  channel 
flow,  have  not  test  filtered  in  this  direction®’^.  We  shall  call  this  form  of  test  filtering 
2D  filtering.  A  typical  test  filter  length  scale  of  twice  the  ‘regular’  LES  filter  length  in 
each  homogenous  direction,  or 

Ax  =  2Ax,  Ay  =  2Ay,  and  Az  =  Az  (17) 


was  used.  In  the  preceding  section,  we  defined  an  explicit  filter  in  the  inhomogeneous 
direction.  This  enables  us  to  test  filter  in  all  three  dimensions,  with  Az  %  2Az  at  all  z 
points.  The  length  scales  Az  and  Az  are  computed  numerically  by  measuring  the  width 
of  a  filtered  delta  function  S{z  —  zjg)  for  both  the  regular  and  test  filters.  The  composite 
length  scale  ratio  A/ A  generally  varies  between  1.85  and  2.15  in  a  typical  application.  We 
shall  call  this  form  of  test  filtering  3D  filtering. 


To  address  the  issue  of  C  and  Pry  spatial  dependence,  we  adopt  the  procedure  of 
Gennano  et  oi.,  and  plane  average  the  tensorial  quantities  in  (16.1)  and  (16.2)  before 
solving 


1  <  LjjMjj  > 

2  <  MijMij  > 


(18.1) 


Pry  =  —C 


<  NjNj  > 

<  KjNj  > 


(18.2) 


This  avera^ng  procedure  is  necessary  to  remove  local  instabilities  that  result  from  point* 
wise  calctilations  of  C .  The  spatial  dependence  of  i/y  and  Pry  is  therefore  restricted  to 
the  normal  directions 


*-r(z)  =  C'(z).A^<15|> 

/cy(2)  =  I^(z)/Pry(z) 


(19.1) 

(19.2) 


with  both  quantities  also  being  functions  of  time.  This  procedure  is  used  for  both  ID  and 
ZD  filtering. 

The  spatial  dependence  of  C  and  Prj*  is  the  subject  of  ongoing  research.  Ideally,  C 
and  Prj<  could  be  functions  of  the  procedure  of  plane  averaging  to  remove 

(x,y)  dependence  has  been  found  to  produce  acceptable  results.  This  is  an  area  needing 
considerable  attention  for  flows  which  are  not  homogeneous  in  these  planes.  Local  aver¬ 
aging  techniques  in  both  space  and  time  are  currently  under  investigation  to  remove  this 
restriction^. 

LES  of  the  Neutral  Ekman  Layer 

The  dynamic  model  was  applied  to  LES  of  a  neutral  Ekman  layer  to  examine  the 
differences  between  2D  and  ZD  filtering.  The  neutral  Ekman  layer  is  that  defined  by 
solving  (3),  (4)  with  gravity  5  =  0.  The  case  considered  here  is  defined  by  the  following 
parameters 

=  B  =  (win|)‘/=,  G  =  (G„o,o),  n  =  (o,o,n,)  (20) 

Coleman  performed  a  direct  numerical  simulation  (DNS)  for  Re  —  500  using  a  128x128x50 
grid  with  length  2S  in  each  of  the  (x,  y,  z)  directions,  where  S  is  the  boimdary  layer  thick¬ 
ness.  In  this  study,  two  LES  grids  with  the  same  lengths  as  the  DNS  grid  were  used, 
a  coarse  32x32x21  grid  and  a  fine  32x32x30  grid.  The  runs  are  labeled  with  a  grid  res¬ 
olution  prefix  and  test  filtering  siiffix,  e.g.  Coarse2D.  Flow  statistics  are  generated  by 
time  averaging  until  the  mean  profiles  satisfy  the  Reynolds-averaged  momentum  equations 
to  within  two  percent  of  the  equilibrium  value.  AH  results  presented  in  this  section  are 
time-averaged. 

The  fidction  velocity,  u«,  and  surface  shear  stress  angle,  /?,  are  the  primary  global 
quantities  used  for  comparisons.  Table  I  gives  the  values  of  u«,  /?,  and  the  amount  of 
resolved  turbulent  kinetic  energy  (TKE)  for  all  LES  runs,  along  with  the  results  from 
Coleman’s  DNS  (since  u*  and  0  are  derived  &om  mean  quantities,  filtering  the  DNS 
results  will  have  a  negligible  effect  on  these  values).  The  results  show  that  3D  filtering 
significantly  improves  the  coarse  grid  calculation.  Indeed,  the  resTilts  from  the  CoarseZD 
calculation  are  better  than  those  from  the  Fine2D  calculation.  The  acctiracy  of  the  LES 
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correlates  well  with  the  fraction  of  TKE  that  is  resolved,  and  in  aU  cases  the  SGS  terms 
are  required  to  model  a  significant  portion  of  the  flow  field. 

The  eddy  viscosity  v^{z)  is  shown  in  Figure  la  for  both  the  near  wall  and  faj  field 
regions.  The  most  striking  difference  is  the  large  disparity  in  uj-  between  ZD  and  2D 
filtering  runs  in  the  outer  regions  of  the  flow,  with  significant  deviation  first  occurring 
aroimd  z"*"  =  30.  Large  values  of  i/y  reduce  the  TKE,  leading  to  a  less  turbulent  flow  as 
measured  by  u*  amd  (laminar  Ekman  layers  have  (3  —  45**,  while  an  ‘infinitely  turbulent’ 
Ekman  layer  would  have  p  —*  O^*). 

Varying  the  grid  resolution  affects  the  profile  of  i/y.  In  Figure  lb,  the  coefficient  C{z) 
in  (19.1)  is  shown.  In  the  near  wall  region  0  <  -■*"  <  30,  C{z)  is  essentially  independent 
of  grid  resolution  and  the  filtering  technique.  This  produces  lower  for  the  Fine2D  and 
FineZD  runs  in  the  near  wall  region,  5. nee  the  value  of  A  is  smaller  on  the  fine  grid  than 
the  coarse  grid.  In  the  outer  region,  however,  the  Fine2D  is  neairly  the  same  as  that  of 
Coarse2Dy  so  that  in  Fine2D  and  Coarse2D  the  same  amount  of  SGS  energy  is  modelled. 
The  same  behavior  is  observed  in  FineZD  and  CoarseZD.  This  indicates  that  the  extra 
resolution  in  the  outer  region  for  fine  grid  calculations  does  not  provide  any  benefit,  and 
that  the  accuracy  of  the  calcxilation  in  the  outer  region  cannot  be  improved  unless  a  finer 
grid  in  the  streamwise  and  spanwise  directions  is  introduced.  Note  that  the  values  of  C{z) 
for  Fine2D  and  FineZD  are  larger  than  their  coarse  grid  counterparts  in  the  outer  region. 

From  Table  I  we  note  that  while  the  FineZD  case  is  the  most  accurate,  it  does  not 
represent  a  substantial  improvement  over  CoarseZD.  We  can  infer  this  from  Figure  la, 
since  the  CoarseZD  and  FineZD  values  of  are  very  close  in  the  entire  flow,  except  for 
the  near  wall  region. 

LES  of  the  Stratified  Ekman  Layer 

In  this  section,  we  shall  examine  the  behavior  of  the  dynamic  subgrid  scale  model  when 
a  mean  temperature  profile  is  imposed  upon  the  Ekman  layer.  The  equations  governing 
the  flow  are  (3)  —  (5),  with  the  gravity  force  g  nonzero. 

The  temperature  profile  will  be  suddenly  imposed  on  the  Fine2D  and  FineZD  neutral 
cases  of  the  preceding  section  at  time  =  0.  The  profile  is  the  S3rmmetric  error-function 
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used  by  Coleman  as  an  idealization  of  the  Convective  Boundary  Layer  (CBL),  and  is  shown 
in  Figure  2.  The  profile  has  both  an  unstable  region,  caused  by  a  surface  heat  flux,  and 
a  stable  region,  due  to  a  capping  inversion  some  distance  above  the  surface.  An  initial 
stirface  lapse  rate  is  selected  and  incorporated  into  a  surface  Richardson  number 


Jlia.i  = 


sro,iD2 


(21) 


Too|G|2 

where  0  subscripts  denote  surface  values.  The  values  of  the  shape  factors  and  surface 
Richardson  number  were  selected  as  =  4JD,  5,  =  5D  «  0.4^,  and  i2io,i  =  —0.05  to 
produce  flow  conditions  consistent  with  the  ‘moderate  heating’  DNS  case  of  Coleman.  Un¬ 
fortunately,  direct  comparison  between  DNS  and  LES  is  not  possible  for  the  stratified  flow 
cases  due  to  a  slight  difference  in  Reynolds  numbers.  However,  correct  interpretation  of 
^he  LES  data  should  be  possible  considering  the  well  known  ‘mixed-layer’  behavior  of  the 
CBL^.  In  addition  to  the  two  LES  involving  the  dynamic  model,  a  third  case,  denoted 
Constant,  was  run  with  constant  values  of  the  Smagorinsky  coefficient  C  —  0.008  and 
Pry  =  0.5  to  illustrate  the  shortcomings  of  such  an  approach  in  the  stratified  Ekman 
layer.  The  definition  of  C  and  the  calculation  of  length  scales  is  somewhat  different  be¬ 
tween  our  LES  and  ‘standard’  finite-difference  techniques  employed  by  the  meteorological 
community^;  OUT  C  corresponds  to  a  Smagorinsky  constant  Cs  »  0.2. 

Figtire  3  shows  the  time  histories  of  global  TKE  for  the  Fint2D,  FintZD,  and  C onsiant 
cases  with  imposed  temperature  profiles.  Time  is  normalized  by  the  global  rotation  time 
scale,  f~^  =  (211^)“^.  We  shall  refer  to  the  first  two  cases  sis  H2D  and  HZD,  respectively, 
since  grid  resolution  is  not  varied.  All  cases  display  an  initial  transient  of  duration  tf  »  0.4. 
HZD  and  H2D  then  approach  a  state  of  quasi-equilibrium,  while  Constant  displays  a 
monotonic  decrease  in  TKE.  Quasi-equilibrium  is  a  well  known  characteristic  of  the  CBL'’^ 
and  is  indicative  of  th'i  rapid  response  of  the  turbulence  on  the  time  scale  reqiiired  to  alter 
the  mean  9  profile.  HZD  resolves  more  TKE  than  H2D,  a  characteristic  observed  in  the 
neutral  Ekman  layer  and  attributed  to  a  smaller  value  of  i^(z).  Time  averaging  is  done 
for  the  period  1.5  ^  t/  <  2.0  to  generate  statistics;  the  slight  unsteadiness  in  the  flow  field, 
along  with  the  small  sample  time,  produces  some  statistical  uncertainty.  Unless  otherwise 
noted,  the  quantities  for  HZD  and  H2D  are  assumed  time-averaged  over  this  interval. 
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Time  averaging  is  not  done  for  Constant,  due  to  the  absence  of  a  quasi- equilibrium  state; 
instead,  instantaneous  profiles  at  tf  =  1.62  are  used. 

Figure  4  shows  the  mean  9  profile  along  with  the  initial  temperature  profile.  Both  HZD 
and  H2D  are  close,  and  both  deviate  substantially  from  Constant.  The  tendency  of  the 
Smagotinsky  model  to  smear  the  capping  inversion  was  first  observed  by  Deardorff^,  and 
is  also  observed  here.  Figure  5  shows  the  buoyant  heat  flux  <  w'd'  >  the  minimum  of 
which  defines  zj,  the  inversion  height  and  a  relevant  length  scale  for  the  CBL.  HZD  and 
H2D  agree  on  the  height  of  both  zj  and  the  location  of  maximum  heat  flux.  The  profiles 
of  9  and  <  w’9^  >  agree  closely  with  accepted  behavior  of  the  CBL,  amd  demonstrate 
that  within  and  above  the  stable  inversion  layer  turbulent  activity  is  largely  suppressed. 
Coleman  calculated  the  convection  velocity,  first  used  by  Deardorff,  as 

w.  =  (gQo^i/Too)^^^  (22) 

and  showed  that,  when  scaled  with  and  zj,  his  DNS  restilts  are  in  excellent  agreement 
with  high-Ae  LES,  laboratory,  and  field  data.  We  will  not  repeat  that  exercise  here; 
instead,  a  comparison  of  global  quantities  for  similarly  heated  cases  is  shown  in  Table  II. 
Good  quantitative  agreement  between  DNS  and  LES  is  observed.  The  quantity  gzj/Too 
is  essentially  constant  in  both  the  DNS  and  LES;  u*  does  not  vary  much  with  Re,  so  we 
would  expect  to  be  lower  for  the  LES  since  the  heating  is  slightly  stronger.  From  the 

calculation  of  w«,  we  note  that  the  time  avera^ng  period  for  the  LES  statistics  corresponds 
to  roughly  2  ‘eddy  turnover’  times  as  defined  by  zj/w*. 

Figures  6a  and  6b  show  vj'iz)  and  C{z)  near  the  wall  and  over  the  bulk  of  the  CBL. 
There  is  a  substantial  difference  between  the  H2D  and  HZD  cases,  with  the  3D  filtering 
procedi:re  producing  much  lower  values  of  and  C,  as  was  the  case  in  the  neutral  layer. 
The  values  of  diverge  much  closer  to  the  wall  than  in  the  neutral  case,  with  significant 
differences  occtirring  at  about  z'^  =  10. 

Figures  7a  and  7b  show  profiles  of  kt(z)  and  the  turbulent  SGS  Prandtl  number,  defined 
as  Prj'  =  i/y/zcy'.  Absolute  values  of  Pry  greater  than  three  are  observed  in  both  HZD 
and  H2D  but  are  not  shown;  these  large  values  occur  when  fcy  fa  0.  For  H2D,  Pr^  attains 
very  large  values  at  points  within  the  CBL  due  the  small  values  of  /cy.  For  HZD,  Pr^  is 
close  to  unity  for  much  of  the  region  0  <  z/zj  <  0.5,  after  which  it  drops  to  Pry  =  0.5  in 
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the  upper  CBL,  and  then  increases  to  Pr’p  »  2.5  in  the  inversion  above  z  =  zj.  Most  LES 
researchers  use  Prj'  values  between  0.33  and  0.5  for  the  CBL^.  One  explanation  for  the 
behavior  of  Prj>  in  this  study  is  that  at  low  iZe,  a  significant  portion  of  the  CBL  near  the 
surface  is  still  influenced  by  the  surface  shear,  and  shear  layers  require  Pry  close  to  1.0. 
Only  away  from  the  surface,  where  the  motion  produced  by  buoyant  eddies  dominates, 
should  Pry  drop.  In  any  case,  it  is  important  to  note  the  variation  of  Pry  with  space 
and  filtering  technique  (though  we  assume  the  H2D  results  to  be  somewhat  erroneous,  as 
will  be  discussed  later).  Large  variations  of  Pry  with  flow  type,  initial  conditions,  and 
even  grid  resolution  have  been  reported  by  Moin  et  al.^,  and  cast  doubt  on  the  practice  of 
assuming  constant  Pry. 

To  the  author’s  knowledge,  the  dynamic  subgrid  scale  model  is  the  first  to  permit 
accurate  LES  of  neutral,  stably  stratified,  and  unstably  stratified  regimes  of  Ekman  layer 
flow  without  the  need  for  modifications  like  (e)  and  (f)  presented  in  the  introduction. 

Discussion 

1.  Effects  of  Inhomogeneous  Filtering  on  the  Dynamic  Model 

As  we  witnessed  in  both  neutral  and  stratified  Ekman  layer  flow,  there  are  marked 
differences  in  the  subgtid  scale  stress  and  heat  flux  depending  on  whether  2D  or  3D 
filtering  is  used  in  dynamic  calculation.  A  sxispected  criterion  for  this  is  as  follows. 
When  the  veloci'  ariance  in  the  inhomogeneous  direction,  w',  contributes  significantly 
to  the  TKE,  3D  filtering  should  be  employed.  This  is  seen  in  Figure  8a,  which  plots  the 
three  components  of  velocity  variance  in  the  neutral  Ekman  layer  LES  FineSD.  In  the 
outer  region,  u'  fts  v  w’  and  the  turbtilence  is  nearly  isotropic.  In  this  ae  region  the 
largest  difference  between  2D  and  ZD  filtering  is  observed. 

If  LES  were  done  on  a  homogeneous,  isotropic  flow  with  equispaced  grid  points  in  aU 
directions,  it  wotdd  be  difficult  to  defend  the  practice  of  abandoning  test  filtering  in  one 
arbitrary  direction  while  retaining  it  in  the  other  two.  Although  the  Ekman  layer  has  an 
inner  region  which  is  quite  inhomogeneous,  its  outer  region,  as  well  as  c  ater  regions  of 
other  external  boundary  layers,  is  nearly  homogeneous  and  isotropic  so  3D  filtering  should 
be  required. 
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The  JETSD  and  H2D  cases  show  a  disparity  in  v^{z)  much  closer  to  the  surface.  This  is 
expected,  since  surface  heating  produces  a  well  mixed,  relatively  homogeneous  flow  with 
values  of  that  can  exceed  v,.  Figure  8b  shows  velocity  variances  for  HZD.  In  the 
neutral  layer,  at  =  30  a  large  disparity  in  C(z)  was  observed;  comparing  magnitudes 
of  velocity  variances  for  FineZD  and  HZD  show  that  in  HZD,  the  magnitudes  aire  much 
closer  at  =  30  than  in  FineZD.  We  would  therefore  expect  a  disparity  to  occur  below 
this  point  for  HZD. 

The  HZD  cadculation  of  the  CBL  is  significantly  more  accurate  than  the  H2D  simula¬ 
tion,  demonstrating  that  2D  filtering  should  be  discouraged  for  the  CBL.  Indeed,  many 
studies  of  the  CBL  using  LES  employ  equispaced  grids  in  all  three  directions  to  tadce  ad- 
vamtage  of  the  ‘quasi-homogeneous’  behavior.  In  such  simulations,  it  should  be  an  easy 
matter  to  test  filter  in  all  three  dimensions. 

A  logical  question  is  ‘Why  has  the  djmamic  model  had  such  success  in  turbulent  channel 
flow  without  ZD  filtering?’  One  possible  explanation  is  that  channel  flow  is  internal;  that 
is,  it  has  a  statistically  steady  state  in  which  the  pressure  gradient  balances  the  total  shear 
stress  gradient  (sum  of  resolved  and  modelled  components).  The  guarantee  of  a  correct 
shear  stress  profile  may  make  up  for  inadequacies  present  in  the  SGS  modelling  technique^. 
In  an  external  boimdary  layer  this  is  not  possible.  . 

U.  Efiects  of  Buoyant  Forcing  on  the  Dynamic  Model 

A  single  subgrid  model  applicable  to  neutral,  unstable,  and  stable  regimes  of  turbulent 
flow  has  been  previously  unrealized.  As  mentioned  in  the  first  section,  several  researchers 
proposed  methods  of  alleviating  the  problem,  which  may  be  characterized  as  a  need  for 
the  eddy  viscosity  to  vary  as  the  stratification  changes. 

The  d]mamic  model  shows  the  ability  to  extinguish  the  eddy  viscosity  in  re^ons  of 
strong  stable  stratification  without  resorting  to  ad  hoc  modifications.  Also  unnecessary  is 
the  additional  expense  of  a  SGS  TKE  equation  to  reduce  the  eddy  viscosity  in  regions  of 
damped  turbulent  activity  (the  solution  of  such  an  equation  involves  further  terms  that 
need  to  be  parameterized).  The  dynamic  calculation  of  ifj>  and  K<p  produces  Prj'  values 
that  are  in  good  agreement  with  accepted  values  in  the  shear  and  convectively  dominated 
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regions  of  the  CBL.  The  failure  of  the  Smagorinsky  model  with  a  constant  C  is  evident  in 
Figures  4  and  5.  In  addition  to  smearing  the  inversion  layer,  Figure  4  illustrates  how  the 
surface  flux  has  been  eradicated.  This  produces  a  poor  heat  flux  profile  that  is  inconsistent 
with  well-known  ‘mixed-layer’  behavior  of  the  CBL.  Conceivably,  some  of  these  deficiencies 
could  be  alleviated  with  the  introduction  of  commonly  used  damping  factors  near  the  wall. 
In  the  spirit  of  model  comparison,  however,  we  have  avoided  ad  hoc  formulations. 

The  widespread  success  of  the  dynamic  model  is  unparalled  and  leads  one  to  believe  that 
the  model  can  be  used  for  other  scalar  quantities.  A  specific  scalar  of  interest  is  moisture, 
with  the  associated  efiects  of  condensation  and  evaporation.  The  introduction  of  such  a 
scalar  into  the  stratified  Ekman  layer  idealization  of  the  Planetary  Boundary  Layer  wiU 
be  the  focus  of  futxire  work. 
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TABLE  I 


Global  Resxilts  for  Neutral  Ekman  Layer  with  Re  =  500 

Grid  and  Filter  u«  x  10^  0  Resolved  TKE 

(DNS  =  1) 

128x128x50 

DNS  6.27  25.4  1.0 

32x32x21 


Coarse  2D 

5.83 

29.7 

0.49 

Coarse  ZD 

6.00 

27.3 

0.64 

32x32x30 

Fine  2D 

5.93 

28.9 

0.58 

Fine  ZD 


6.02 


27.3 


0.70 


TABLE  II 


Global  Results  for  Stratified  Ekman  Layer 


Grid  iEtO.i 

128x128x50 
DNS  -0.04 

32x32x30 

EZD  -0.05 


if  u*  X  10^ 


2.1  7.48 


1.5  -  2.0 


/?  (deg) 

28.1 


0.93 


6.93 


29.2 


0.95 


Figure  la.  Eddy  viscosity  for  neutral  Ekman  layer. - ,  FineZD^ - ,  Fine2D 

CoerteZD, - ,  CoarseZD 
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Figuxe  2.  Initial  temperature  profile. 


Figure  3.  Time  history  of  globally  integrated  TK£  for  stratified  Ekman  layer. 
- .  H2D. .  .  CoMtant 


Figure  5.  Mean 
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Figuie  6b.  Smagoiinskj  Coefficient  for  stratified  Ekman  layer. 
.  ,  Constant, 


,  HZD, - ,  H2D, 


Figure  8a.  Velocity  variance  for  neutral  Ekman  layer  ca»e  FineZD 
(streamwiae), .  ,  v  (spanwise), - ,  w  (vertical). 


Figure  8b.  Velocity  wiaace  for  stratified  Ekman  layer  case  HZD 
(streamwise), .  ,  v  (spanwise), - ,  w  (vertical). 


